Light Quark Masses to All Order of Chiral Expansion 
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I. INTRODUCTION 



The quark masses are some of basic parameters of the standard model. There are various hadronic phenomenologies 
relating to the hght quark(u, d, s) masses. For instance, they break the chiral symmetry of QCD expUcitly. m„ — 
breaks isospin symmetry or charge symmetry, and (rn„ + TO(j)/2 — TOs breaks SU(3) symmetry in hadron physics 
, respectively. However, in QCD the masses of the light quarks are not directly measureable in inertial experiments, 
but enter the theory only indrectly as parameters in the fundamental lagrangian. The purpose of this paper is to 
[ study light quark masses at energy scale /i — rup in a non-perturbative way. 

I In general, at low energy information about the light quark mass ratios is extracted in order by order by a rigorous, 
' semiphenomenological method, chiral perturbative theory(ChPT) To the first order, the results, mg/md = 19 

. and rriu/niii = 0.556, are well-known Many authors have studied the mass ratios up to next to leading order 
' of the chiral expansion. Gasser and Leutwyler first obtained ms/md = 20.2 and mu/md = 0.554 Then Kaplan 
Oh, and Manohar extracted nis/md =15 to 23 and rriu/md =0 to 0.8 with very larger error bar These values have 
f-^ • been improved to ms/md = (20.5 ± 2.5), mu/md = 0.52 ± 0.13 by Leutwyler fli^^o ms/md = 18, mu/md — 0.66 
D ', by Gerard j?), and to ms/md = 21, mu/md = 0.30 ± 0.07 by Donoghue et. respectively. Finally, Leutwyler 

> analysis previous results and obtained ms/md — 18.9 ± 0.8, mu/md = 0.553 ± 0.043 Q. So far, however, the study 
^ on light quark masses in framework of ChPT is limited by the following shortages: 1) In this framework, to obtain 
light quark mass ratio beyond the next to leading-order is very diffcult, since more and more free parameters are 
included with raising of perturbative order. 2) Due to Kaplan-Manohar ambiguity the defination of the light 
quark masses in ChPT is not uniquely. The reason has been pointed out in ref. ||8|] that, at the next to leading-order, 
QCD renormalization is mass-dependent, and the symmetry can not distinguish renormalized quark masses from 
"bare" quark masses. 3) In framework of ChPT, we can only obtain light quark masse ratios. For obtaining the 
individual quark masses, other approachs, such as QCD sum rules or lattice calculation [|l2| are needed. These 

shortages are urgently wanted to be improved by theorectical studies of QCD and experiment. 

In ref. we have the constructed chiral constituent quark model(ChCQM) including the lowest vector meson 
resonances following the spirit of Manohar-Georgi model jl^ . This model provides a formulation to perform rigorous 
field theory calculation at energy scale lying between ChPT(/i ^ 0.5GeV) and chiral symmetry spontaneously bro- 
ken(CSSB) scale(^ ^ 1.2GeV), and a successful description on physics in this energy region Jl3| , p^ . Up to 0(p*), low 
energy limit of the model agree with ChPT well. Thus this model can be treated as an approach to extend ChPT 
investigation inspired by QCD. The most important advantage of this approach is that we can perform calculation 
beyond the low energy expansion in ChCQM, and only fewer free parameters are required. In general, there are three 
types of expansion working at low energy. They are momentum expansion, light quark mass expansion and 
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expansion ]16| |. In ref. |13|Jl5| |, we have provided a rigorous method to perform calculation to include all order terms 
of the momentum expansion and up to the next to leading order of expansion. In this the present paper, we will 
extend this method to reflect all order information of light quark mass expansion. 

The Kaplan-Manohar ambiguity of ChPT has caused many debates. In particular, due to this ambiguity, the 
authors of ref argued that the observed mass spectrum is consistent with a broad range of quark mass ratios, 
which specially includes the possibility niu = 0. However, it is disagreed by other authors [l^J^. This problem can 
also be discussed in ChCQM. Since ChCQM is an effective approach with features of low energy QCD, light current 
quark masses are defined uniquely in ChCQM, that are just renormalized "physical" masses of u, d, and s quarks. In 
principle, therefore, there is no Kaplan-Manohar ambiguity in ChCQM. Light current quark masses can be determined 
uniquely via meson spectrum. 

In ref. @, we studied the isospin breaking process lu tt+tt and determined isospin breaking parameter — 
m„ = (3.9 ± 0.22)MeV at vector meson energy scale. This result together with meson spectrum provide so much 
information that we can obtain not only light quark mass ratios but also individual quark masses. 

In ref. [T^ , we have shown that the chiral expansion at vector meson energy scale converge slowly. In particular, 
we have also pointed out that, if we neglect strange quark masses, the chiral expansion at (/)(1020) energy scale will 
be divergent! It implies that nis play very important role at (/)-physics. In ref. jl3Ulq| we have successfully studies the 
chiral expansion at nip and energy scale. In order to extend this study to if* (892) and 0(1020), individual quark 
masses are neccessary. It has been recognized that the light quark masses obtained in different approaches are with 
larger difference. Thus we have to extract quark masses by this formalism itself. It is one of our goals. 

In general, in ChCQM the light quark masses can be extracted not only by pseudoscalar meson spectrums, but also 
by the lowest vector meson resonance spectrums. However, one-loop effects of mesons will contribute to vector meson 
masses, and calculation on one-loop effects of mesons is related to the chiral expansion at vector meson energy scale. 
As shown in ref. [T^ ] , this relation is very complicate. It makes that the relationship between vector meson spectrums 
and the light quark masses are also very complicate and indirect. Thus in this paper we will extract information 
about the light current quark masses from pseudoscalar meson spectrums and their decay constants. The vector 
meson spectrums will be predicted by this formalism in other paper. 

The contents of the paper are organized as follows. In sect. 2 we review the basic notations of the chiral constituent 
quark model with the lowest vector meson resonances. In sect. 3, masses and decay constants of pseudoscalar meson 
octet are calculated. The results will including all order information of the light quark masses. In sect. 4, one-loop 
effects of pseudoscalar mesons and renormalization are discussed. The numerical results are given in sect. 5 and a 
brief summary is included in sect. 6. 



II. CHIRAL CONSTITUENT QUARK MODEL 

The s imp lest version of chiral quark model which was originated by Weinberg |l[] , and developed by Manohar and 
Georgi p4[ provides a QCD-inspired description on the simple constituent quark model. In view of this model, in the 
energy region between the CSSB scale and the confinement scale (Aqcd ^ 0.1 — 0.3Gey), the dynamical field degrees 
of freedom are constituent quarks(quasi-particle of quarks), gluons and Goldstone bosons associated with CSSB (these 
Goldstone bosons correspond to lowest pseudoscalar octet). In this quasiparticle description, the effective coupling 
between gluon and quarks is small and the important interaction is the coupling between quarks and Goldstone 
bosons. In I we have further included the lowest vector meson resonances into this formalism. At chiral limit, this 
model is parameterized by the following chiral constituent quark lagrangian 

= iqi/d + /T + 5^^75 - i/V)q - mqq + V^UV^U^ > +Knl < F^V^ > . (1) 

Here < ... > denotes trace in SU(3) flavour space, q = (qu,qd,qs) are constituent quark fields. denotes vector 
meson octet and singlet. Since in this paper we only focus on pseudoscalar meson spectrums and decay constants, we 
will neglect vector meson fields in the following. The and are defined as follows, 

= liCKd^ - ir^)^ + - zl^)^}, (2) 
and covariant derivative are defined as follows 

V^[/t = 0^[/t - tlpU^ + tU^r^ = -2e^A^e^ (3) 
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where Ip, = v^^ + a^^ and = Vp, — Ufj, are linear combinations of external vector field and axial-vector field a„, £ 
associates with non-linear realization of spontanoeusly broken global chiral symmetry introduced by Weinberg [[l9[ . 
This realization is obtained by specifying the action of global chiral group G = SU{3)l x 517(3)^ on element ^(<i>) of 
the coset space G/SU{3)^: 

m ^ 9Bm)hH^) = h{<^>)a<^)g[, 9L,9ReG, h{^)eH = SU{3),. (4) 

Explicit form of ^($) is usual taken 

^($) =exp{iA°$°(a;)/2}, ^2 (5) 

where the Goldstone boson are treated as pseudoscalar meson octet. In rcf. we have shown that the la- 
grangian( |l|) is invariant under Ggiobai x Giocai- 

Distingushing from some versions of chiral quark models, there is a kinetic term of pseudoscalar mesons in la- 
grangian( 0). Therefore, the kinetic term of pseudoscalar mesons generated by one- loop effects of constituent quarks 
can be renormalized. Note that there is no mass term of pseudoscalar mesons in eq. (^. All of these are due to 
the basic assumption of the model In the other words, in this energy region, the dynamical field degrees of 

freedom are constituent quarks and massless Goldstone bosons(pseudoscalar octet) associated with CSSB. Masses of 
pseudoscalar mesons will be genarated by quark loops as current quark mass parameters emerge in the dynamics (to 
see below). In ref. jl3| we have fitted the parameter qa = 0.75 via /3-decay of neutron, and m = 480MeV via low 
energy limit of the model. It has been also pointed out that the value of gA has included effects of intermediate 
axial- vector meson resonances exchanges at low energy. 

The light current quark mass-dependent term has been introduced in ref. ]l5[ | based on requirement of the chiral 
symmetry, 

- Iqi^^xe + ^X^Oq - f - ex^0759, (6) 

where x = s + ip, s = Scxt + M, A4 — diag{m„, md, TOs} is light current quark mass matrix, Soxt and p are scalar 
and pseudoscalar external field respectively. Eq. (^) will return to standard quark mass term of QCD, '!/'A^'0, in 
absence of pseudoscalar mesons at high energy for arbitrary k. Therefore, although the light current quark masses 
are defined uniquely in this formalism, the symmetry and the constrains of underlying QCD still can not fixed the 
coupling between pseudoscalar mesons and constituent quarks, k will be treated as an initial parameter of the model 
and be fitted phenomenologically. 

To conclude this section, the ChCQM lagrangian with light current quark masses is 

= iq{/d + /r + 3^,^75)9 - mqq - ^qi^^M^ + CX^O^ " f ^(C^X?^ " ^X^Ol^Q 

+ ^ < V^C/V^C/t >, (7) 

where vector meson fields have been omitted. The effects of isospin breaking due to inequality of light quark masses 
will be generated via constituent quark loops. For example, it generated masses of pseudoscalar mesons and splits 
their decay constants. 



III. QUARK LOOPS 



In this section, we will calculate pseudoscalar meson masses and decay constants induced by one-loop effects of 
constituent quarks. They are the leading order in expansion. 

In this framework, the effective action describing meson interaction can be obtained via integrating over degrees of 
freedom of fermions 

e'^^-ff = J VqVqe I '^""'^''''''^ =< vac, out\in,vac >Ext, (8) 

where < vac, out\in, vac >Ext is vacuum expectation value in presence external sources. The above path integral can 
be performed formally, and many methods, such as heat kernel manner pO| , pl| , have been used to regulate the bilinear 
operator yielded via path integral. By using those formal integral methods, however, the explicit calculations of high 
order contributions to the chiral expansion are extremely tedious. In ref. [p^ we have provided a convenient method 
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to evaluate effective action via calculating one-loop diagrams of constituent quarks directly. This method can capture 
all high order contributions of the chiral expansion. 

In interaction picture, the equation( ph is rewritten as follow 



ni(0) + E; / ^•■•7^n„b,,.,p„_,), (9) 



n=1 



where Tq is time-order product of constituent quark fields, £^ is quark-meson interaction part of lagrangian( |^), 

nn(pi, • ■ • ,Pn) is one-loop effects of constituent quarks with n external sources, pi,P2i ■ ■ ■ ,Pn are four-momentas of n 
external sources respectively and 

n„(pi, • • • ,p„_i) = J d^PnfiniPir ■ ■ ,Pn)S'^ipi - P2 Pn) ■ (10) 

To get rid of all disconnected diagrams, we have 

oo 
n=l 

^i = ni(o), (11) 

f d-^pi d^pn^i / ^ o^ 

= y (2^ • • • 72^""^^^' • • • - 

Hereafter we will call Sn as n-point effective action. 

For the purpose of this paper, the lagrangian( ^ can be rewritten as follow 

C^{x) = Cgix)+C^^\x)+C[Tr{x),Mx)]+C[K^{x)]+C[K''{x)l (12) 

where Cq and £2*^'' are free field lagrangian of constituent quarks of pseudoscalar meson respectively, 

Cq = qi{i/d - mi)qi, fhi = m + mi, 

i—u,d,s 

= —d^.^'^d^"^" + Aald^'^'' + ■■■, a = 1,2,- ••,8, (13) 



where 



A"<J>'^(x) = ^/2| TT- -^ + ^ if" 

if- -^,vs 



A' 




£[<i>(x)] denotes quark-meson coupling lagrangian 

9A r/c, _+ , n^+\.-...u.. J , . 9A 



L\k, rys] = ^[(5^7^+ -t- 2a+)i27'^75d c.c] + ^ E + 2A«)*7^759, 

i—u.d,s 



%—u,d,s 
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+ + md)7r+7r (uu + dd) + - ^ miP^qiqi, (14) 



2 

i—u.d,s 



C[K^] = + 2A+)S7''75S + c.c] + + m,)(i^+u75S + c.c.) 

+ -(m„ + ms)K^ {uu + ss), 



2 

i - - 1 

H — ^K,{md + ms)iK°d^5S + K°d-y5s) + - {ma + ms)K°K"{dd + ss), 
V 2 2 

where c.c. denotes charge conjugate term of pervious term, 

1 12 

Pu^TT3 + —1]s, Pd = -713 + —^VS, ^ ^"Tf'^*' ^^^"^ 

and 4'^; etc. are axial- vector external fields correponding these pseudoscalar meson fields. From eq.( |lj) 

we can see that the JsT^-quark coupling and ^fC^-quark coupling are similar to 7r^-quark coupling. Thus we only need 
to calculate masses and decay constants of tt^, 7r° and r]^. Then masses and decay constants of can be obtained 
via replacing by in one of tt^ , and masses and decay constants of can be obtained via replacing to„ by nig 
in one of tt^. 

The one-point effective action is generated by tadpole-loop of constituent quarks. Calculation about this tadpole- 
loop contribution is simple. 

iSi[n,r]s] ~ 2 J ^'*'''{(™" + TOd)7r^(.T)7r~(x) < 0\T{u{x)u{x) + d{x)d{x))\0 > 
+ m^Plix) < 0\Tq,{x)q^{x)\0 > 

i—u,d,s 

= -\j d^x J ^{[{m^ + md)7:+{x)7r'{x)+m^Pf{x)]Tr[S^^\k)] 

+ [{mu + md)n+{x)n-{x)+mdP^{x)]Tr[Sf{k)]+mgP^{x)Tr[SP{k)]}, (16) 

where Tr denotes trace taking over color and Lorentz space, 5^^''' is propagator of constituent quark fields, Sp^ (k) — i{/ 
k — fhq + ie)~^. In terms of dimensional regularization, we can integrate over internal line momenta k in the above 
equation. The result is 



z5i[7r,ry8] = -^^p^r(l-|)z J d^x{{^Y^^rhl[im^ + md)7r+n- + m^P^] 



,2 ,,2 



+ (^2 )'/'m3[(m„ + mrf)7r+^- + m,P|] + (-^)^/2^3m,p2}. (17) 



mi 

Defining a constant Bq to absorbe the quadratic divergence from loop integral, 

?2 AT ,,2 



§s.-^(£.)"^r(i (IS, 



we have 



S'i[7r,778] = J d'^xCi[Tr{x),ris{x)] 

Ci[n,m]^-^Bo{{xl + xl){m^+md)Tr+Tr- + ^ xfm.P^^} 

i—u,d,s 



■m^{{xl\iixl+x^liixl){mu + md)n+n + ^ m^xf Inxfi^^} (19) 



i—u,d.s 
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where Xq = mq/m. 

The two-point effective action concerning to masses and decay constants of charge pion can be obtained as follow, 



+ i(m^ - TO^) - t(l - t)q^ 



{An) 



+^r(2-|)A^5^(m.+mrf)* 1 d^^, j _0_e^«-[(^g^vr+(q) + 2a+(g))a^7r-(x) + c.c] 

xy dt[m„ + i(md-m„)] f ^ " ^^^7 ^ —\ (20) 

x{r(l - + - ml) - t(l - i)q^\ - r(2 - ^)[m„2 _ ^^^^ + t(^2 _ _ 2i(l - 

There are both quadratic divergence and logarithmic divergence in the above effective action. The quadratic divergence 
can be canceled by constant defined in eq.( |8|), and the logarithmic divergence can be conceled via defining 

Fn Nc . n , D.Anu^,, , , 

i-T6 + (4^^^™^(^-2)(^)^- (^1) 



g is an universal coupling constant of this model. In ref. |13|, it has been determined as 5^ = ^ by the first KSRF 
sum rule p2|. 

Then eq.(]20|) together with eq.( |9|) give 0{Nc) effective lagrangian containing the terms linear or quadratic in the 
charge pion fields as follow 

A (22) 



where 



3 

F'^{mu,md) = F^ + -g'^gKniu + md){Am + niy, + m^) + ig"^ ng Ai^u + rnd){rhu + rha) 

--^gAigAifhu + fhd) + 2n{mu + m<j)] j dt[fhu + t{md - m„)] In {xl + t{xl - xD) 

- 3 3 

f{mu,fnd) = F^ + -g'^g\[m^ + md){Am + m„ + m^) + -g'^ngA{mu + md)[mu + m^) 
-^5A[ffA('7!.„ + rhd) + K(m„ + m^)] / dt[TO„ + <(to<j - m.„)] In (x^ + i(a;^ - xD), 
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_ A^^. I' ^^^2^, _ ^^^^ ^ „ rnl)] In (x^ + t{xl - xD), (23) 

Z ZTT _fQ Jg 



a(g ;m„,md) = -^gA[5A(m„ + m^) + K(mu + md)] / dt[m„ + i(md - m„)] In 1 _ j . , , _ 2 -2^/' 
Nc _ _ - ( ^(1 - 



27r2 i^2l y ^2+i(^2_„^2) 



Jo \ ml + t{mj-ml)J 

fhufhd 



It should be pointed out that a{q'^;mu,md) is order q^ at least and (3{q^]mu,'md) is order 5* at least. Since in this 
paper we focus on pseudoscalar meson spectrums, these high order derivative terms should obey motion equation of 
pseudoscalar mesons. In momentum space, the motion equation of physical pseudoscalar mesons is generally written 

{q^ - mlMq) = -lUq'' A^;f\-q), (24) 

where and are physical mass and decay constants of pseudoscalar, e.g., = 135MeV and = 185.2MeV. 
Due to this motion equation, we have 

a{q^ ■,m.u,md)a+{x)-K^ [q) = a{ml;mu,md)a+{x)TT^{q), 

P{q'^;mu,md)TT+{q)TT^ {-q) = /3(m2 ; m„, mrf)7r+(q)7r"(-q) - -g^/^/3'(m^; m„, TO<j)(a+(a;)7r" (q) + c.c), (25) 

where 

dq 

Thus eq.( p^ ) can be rewritten 



I3'{ml;mu,md) = — /3'((?^; m„, TO<j)|q2^„2 . (26) 



+{^^^ + ^ ^^^f^^'^^l,m.,^,)yiatd^^- + C.C.). (27) 

^ z yniii^Tfid) 

The two-point effective action concerning to masses and decay constants of neutral pion and ?78 can be evaluated 
similarly. However the decay constants for neutral mesons cannot be extracted directly from the data. It means that 
the decay constants for neutral mesons cannot be used to determine light current quark masses. Therefore, in this 
paper we do not need to evaluate the decay constants for neutral pion and 773. The effective action containing the 
quadratic terms of neutral pion and r^g is then 

iS'2[7r°,778] 

= f/(0(S^ ^ <l'PA^)PA-<l)Tr\l''T.S'P{k~q)ri,S'P{k)] 

j—u,d,s 

+ '-gAK J2 / St {iq,)PMPj{-q)Tr[ri^S'-f{k - q)l,S^P {k)] 

j—u,d,s 

Y E I (2^)4 ^t^^ PMPA^'imToSPik-qh.sPik)] 



j—u.d,s 



j—u,d,s 
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j—U,d.S \ LL \ 



x{r(l - -)[mj - t{l - t)q'] + 2t{l - t)q'T{2 - -)}. (28) 

The divergences in the above effective action can be canceled by eqs.( |l^) and (21). Then from eqs.( |l9|) and ( [2^ ) 
we can obtain 0{Nc) effective lagrangian describing two-point vertex of tt^ and rjg, as follow 



'i—u,d,s 



where 



,2 



-^c / 2- 2io - I 22\i 2 

-^(Sa'^i ^dAKmiUii + K mj Inx, , 



Mf(TOO = (So + lna;2)TO,a;2(a;, '-k^), (30) 
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iV /■"'^ t(l — t)a 

P.,{q^;mi) = TT^^l^i^m} + (.g^m^ + 2Km,m, + 2K'mj[fh}q-^ - 3t(l - t)]) / rfHn (1 - ^ _ )}■ 

Since auxiliary fields Pi do not lie in physical hadron spectrums, the equation of motion ( p^ ) can not be used in (|2£ 
simply. In section 5, we will use propagator method to deal with the terms with high power momenta in (|2^) and 
diagonalize tt^ — r/s mixing. 



IV. MESON LOOPS AND RENORMALIZTION 



Purpose of this section is to evaluate one-loop effects of pseudoscalar mesons. Due to parity conservation, there are 
only tadpole diagrams of pseudoscalar mesons contributing to masses of decay constants of 0^ mesons(fig. 1). Since 
ChCQM is a nonrenormalizable effective theory, it is very difficlut to calculate meson loop effects completely. However, 
we can expect that, in mass spectrums and decay constants of pseudoscalar mesons, the dominant one-loop effects 
is generated by the lowest order effective lagrangian, because neglect is O(to^) and suppressed by expansion. 
Fortunately, the formalism is renormilzed up to this order. 

The lowest order effective lagrangian is well-known 

16 

p2 op2 

= ^< [A«,A^][A",A^] > +^ < y'X^'i^xH + ^he) >, (31) 

where A°(a — 1, 2, • • • , 8) are Gell-Mann matrices, x ~ 2i?oX and the following SU(N) completeness relation have been 
used 

V < X'^AX'^B > = -— < AB > +2 < Ax B >, 



E < >< > = 2<AB>-'^<A><B>. (32) 

a=l 

To evaluate the one-loop graphs generated by this lagrangian, we consider the quantum fluctuation ^{x) — ip°'{x)X'^ 
around the solution U{x) = f^(x) to the classical equations of motion, 
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U = ^e^^C. (33) 
Substituting cxpansion( |3^ ) into £2 and retaining terms up to and including ip'^ one obtains 

£2 - £2 + ^{d^^'^d^ip'^ - mlcp-^-) '^<[V, Am][^, A^] + l{cp, ipjih^^+^^X^^ - ^M) >, (34) 

where we have omitted some terms which do not contribute to masses and decay constants via one-loop graphs. The 
contribution of tadpole graphs can be calculated easily 



4'''^ = -J<[A^A^][A^A^] + i{A^A«}(^Vc+M-2X)> J ^ 

1 n rn^^ ni^„ 1 



(27r)4 fc2 - ml + ie 



K^- T?^i^^) < [A^ AM][A^ A^] + -{A^ An(eVe + e^xe^ - ^m) >, (35) 



where 



Comparing eq.( |35| ) and eq.( ^ ) we can see that the divergence A can be absorbed by free parameters Fq and Bq. 
Thus the sum of tree graphs and tadpole contribution is 

4 = 41(1-3^^)^^ < [A^A,,][A^A^] > +2^(1 - -M^) < A'^A'^(Cxte + e^xe^) >, (36) 

where 

2 2 

In lagrangian( Bra) we appoint that m^p — m^r when a = 1,2,3, rn,^ = when a — 4,5,6,7 and rriip = m^^ when 



V. LIGHT QUARK MASS DETERMINATION BEYOND THE CHIRAL PERTURBATION EXPANSION 

For extracting (™^+ )qcd from experimental data, the electromagnetic mass splitting of A'-meson is required. The 
prediction of Dashen theorem |2^], [m^^ — ™^o)e.m. = 1.3MeV, has been corrected in serveral recent analysis with 
considering contribution from vector meson exchange. A larger correction is first obtained by Donoghue, Holstein 
and Wyler |2^, who find [m^^ — m^„)e.m. = 2.3MeV. Then Bijnens and Prades ||2^, who evaluated both long- 
distance contribution using ENJL model and short-distance contribution using perturbative QCD and factorization, 
find (to^+ —^j^o )e.m. — 2.4±0.3MeV at ^ — nip. Gao et.al. also gave (jn^^ ~™jro )e.m. — 2.5MeV. Baur and Urech 
psf , however, obtained a smaller correction, (w^^ — ™^o)e.m. — 1.6MeV at /i = nip. In addition, calculation of lattice 
QCD found (to^+ — ni^g)e.m. = l-9MeV. These estimates indicate that the corrections to Dashen theorem are 
indeed substantial. In this the present paper, we average the above results and take ni^^ — ni^g)e.m. = 2.1 ± O.lMeV 
at energy scaleu = nip. 

From eqs. ( H) and (|3|), the masses and decay constants of koan and chrage pion can be obtained via solve the 
following equations 

"^1+ ^ / ^ ° r{M^(m„,md) + /3(m^;m„,TOd)} 

= F2^^° A Ml{niu,ms) + (3iml;mu,ni,)}, 



f2 




md) 






i^|(m„. 


nis) 


Fi 




F^imd, 


TOs) 


flirriu, 


md) 


F^imu, 


rud) 



{Ml{nid,ms) + (i{nil;md,nis)}, (38) 



_ j^y'^u,''''d, , a{nil;niu,md) , F^ ^ r'(^2.^ ^ . 

J" R[mu,md) I' ^(niu,md) it ^{niu,md) 
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4+ = T^T 7 + ^ + 



P2 
^0 



fKP'{'m'1;md,ms), 



Fj^{md,ms) F^{md,ms) 2F^{md,ms) 
where subscript "R" denotes renormalized quantity, 



F^imu, md) 












= F'^{mi,ms) ^ 








Ml(mu,md) 


= M'^{mu,md) 


- So(m„ 










- Bo{mi 4 


3 3 


5 


flimu,md) 


^P{m^,md)~ 








fl{m„ms) 


= P{mi,ms) - 









(i — u, d) 



(39) 



Here the quantity /x^ is defined in ( |37| ) and depend on the renormalization scale /z. It has been recognized that the 
scale plays an important role in low-energy QCD. In this formalism, many parameters, such as light current quark 
masses, constituent quark mass m, axial- vector coupling constant gA, are scale-dependent. The characteristic scale 
of the model is described by the universal coupling constant g which is determined by the first KSRF sum rule at 
energy scale fj, = nip. Hence we take = rup = 770MeV in /i^, and the physical masses and decay constants, however, 
should be independent of the renormalization scale. 

In order to obtain the masses of tt^ and rf , the tts — 773 mixing in eq.( |9|) must be diagonalized. Eq.(p9|) together 
with eq. (E^) lead to the quadratic terms for the ir^ and ryg are of the form 



<S'2[7r3,778] 



A ,1 

(27r)4^2 



(40) 



where 



1 



p2 _ _ 2 



Miiq') 



^0 r)i#2 



3K 



2 {M^(m„) + Mj{md) + 4.M^{ms) + /3„((j^; m„) + (3d{q^;md) + 4/3s((j2; m,) 
2 . . , 3 1 



-2Bo(m + 2m,)(2/z^, + -/i^J - 2Borh{--fi^ + fi^ + -^J.^s)} 



(41) 



1 



-^=^-^{M^(m„) - Mdirrid) + /3«(g^; mu) - I3d{q^]md) - Bo(m„ - md){-fi.^ + Mk- + ^Miys)} 
, 2F2(m„)-F|(m,) 



V3F3F8 
with m — (m„ + md)/2 and 



Fl = F2(™„) + Fjimd) - F2(2m. + Mk): 

Fl = \{Fl(m^) + Fj(md) + ^F^A^s)} - 3^^oV. 



(42) 



Due to 7r3 — ryg mixing, the "physical" propagators of tt'^ and ryg are obtained via the chain approximation in 
momentum space 



q^-m^o+ii- q^-Mi{q^)+ie [q^ - Miiq^) + ie){q'^ - M^q^) + ief 

i 
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g2 - r7i2^ + ze - Mi{q^) + le ^ [q^ - Ml{q^) + ie){q^ - Ml{q^) + itf 

i 



(43) 



Then the masses of tt" and rjs are just solutions of the foUowing equations, 



mlo-Miimlo)' 



mis = Mi {ml ) + J^'i'?};^ , . (44) 



In eqs.( 23) and ( ^0|), the parameters k, Fq and Bq are still not determined. In order to determine them and three 
light quark masses, six inputs are required. In this paper we choose — 185.2 ± O.SMeV, f + — 226.0 ± 2.5MeV, 
m^o = 134.98MeV, m^^^ = 497.67MeVand (to^+)qcd = 491.6±0.1MeV. Ajiother input is mrf-m„ = 3.9±0.22MeV, 



which is extracted from lu tt+tt" decay at energy scale fi = rrip in ref. |15|. Recalling m = 480MeV, gA = 0.75 and 
g = TT^^ for Nc = 3, we can fit light quark masses as in table 1. 

In table 1, the errors of results are from uncertainties in decay constant of . electromagnetic mass splitting of 
if-mesons and isospin violation parameter rud — niu respectively. The first column in table 1 we show the results for 
/ = 223.5MeV. The third column corresponds the center value of / ^ , 226.0MeV and the fifth column corresponds 
= 228.5MeV. Then from table 1 we have 

m, = 160 ± 15.5MeV, nid = 7.9 ± 2.7MeV, m„ = 4.1 ± 1.5MeV, 

TTi TTi 

— = 20.2 ±3.0, — = 0.5 ± 0.09. (45) 

md md 

Here the large errors are from the uncertainty of /^^^ . From table 1 we also have {mT^+ — m„o)QQYi — — 0.25MeV(in this 
paper the contribution from tts — rj' mixing is neglected). This result allows that the electromagnetic mass splitting 
of 4.8MeV. 





Fit 1 


Fit 2 


Fit 3 


Fit 4 


Fit 5 




185.2 


185.2 


185.2 


185.2 


185.2 


■1 K + 


223.5 


224.1 


226.0 


227.8 


228.5 




134.98 


134.98 


134.98 


134.98 


134.98 


„ 


497.67 


497.67 


497.67 


497.67 


497.67 




491.6 ±0.1 


491.6 ±0.1 


491.6 ± 0.1 


491.6 ±0.1 


491.6 ±0.1 


Biuj TVTV) 


1.95% 


1.95% 


(2.11 ±0.20)% 


(2.21 ±0.25)% 


(2.21 ±0.30)% 


rris 


144.2 


150.0 


161.8 


172.0 


175.4 


rrid 


6.26 


6.83 


7.94 ± 0.07 


8.94 ±0.1 


9.53 ±0.11 


niu 


2.56 


3.13 


4.10 T 0.07 


5.03 T 0.1 


5.63 T 0.11 


rriu + md 


8.92 


9.96 


12.04 


13.97 


15.16 


rrid ~ ruu 


3.7 


3.7 


3.84 ±0.14 


3.9 ±0.2 


3.9 ±0.22 




224.0 


224.7 


226.6 


228.4 


229.0 




134.74 


134.74 


134.73 


134.73 


134.73 




608.8 


616.9 


628.8 


636.8 


641.0 


K 


0.5 


0.4 


0.3 


0.25 


0.2 


Fo 


156.95 


156.7 


156.1 


155.6 


155.3 


Bo 


2141.1 


1886.6 


1551.1 


1328.6 


1216.8 



TABLE I. Light current quark masses predicted by the masses and decay constants of 
pseudoscalar mesons at energy scale jj, = rup. Here k is dimensionless, other dimensionful 
quantities are in MeV, and f denotes input. The formula for branching ratio oi lu ^ tt+tt" 
can be found in ref. [15]. It determines the isospin breaking patameter nid — mu- 



ll 



VI. DISCUSSION AND SUMMARY 



In this paper we study information on light quark masses at energy scale fi = nip in the framework of chiral 
constituent quark model. The analysis of quark masses beyond the next to leading order in the chiral expansion is a 
challenging subject. The attempt is stoped in ChPT due to the difficulties mentioned in Introduction. In approach of 
quark model, however, it is not necessary to treat the light current quark masses as the quantities used to construct a 
perturbative expansion. In addtion, since the light current quark masses can be defined uniquely, the Kaplan-Manohar 
ambiguity is avoided in this formalism. From eq. ( ^ ) rriy , 7^ is also confirmed in our results. This conclusion supports 
viewpoint of ref. jlj,^ but disagree with one of ref. ||l7|] . 

Our results are yielded by a non-perturbativc method and contain complete information on light quark masses. Not 
only quark mass ratios but also individual light quark masses are obtained. The results agree with one obtained by 
other approachs(e.g., QCD sum rule or ChPT) well. 

Finally, it is interesting to expand non-perturtation results (^) up to the next to leading order of the chiral 
expansion and compare with one of ChPT. If we neglect the mass difference nid — Wu, up to this order the masses of 
decay constants of koan and pion read 



ml 



B,{rn^ + m,){l + - + [-(3 - ^ + -^(_ - - _ - 

B^im + 7..){1 + + [-(3 - + - ... - f ~ ^^^)]^}' (46) 

U = Fa{l - AtTT - ^fJ-K + ^9a 



ml. 



^ , 3 3 3 3 2 mirh + rris, 
^o{l - ^/"x - ^Mk- - gM^8 + ^9a }■ 



Comparing the above equation with one of ChPT we predict the O(p^) chiral coupling constants, L4, L5 Lq and 
Ls, as follow 



3to 2 



L4 — Le — 0, L5 — 2u 9a 

T _ ^0 (o ,.2n , 3?77 TO g\ ^0 „2 X , ^ 



Numerically, inuptting k. = 0.35±0.15, Fq = 0.156GeV and Bq = 1.6±0.4GeV, we obtain L5 = (1.6±0.3) x IQ-^ and 
Lg = (0.7±0.5) X 10"^. These value well agree with one of ChPT, = (1.4±0.5) x lO'^ and Lg = (0.9±0.3) x lO'^ 
at energy scale /i — nip. This fact can interpret why the light quark mass ratio in eq. (^ ) are close to results extracted 
from ChPT by Leutwyler |^ . The above results indicate that the contribution from scalar meson resonance exchange 
is small. In fact, in hardon spectrum, there is no scalar meson otect or singnlet which belong to composited fields of 
qq. Thus it is a ad hoc assumption to agrue some low energy coupling constants, such as and L^, receiving large 
contribution from scalar meson exchange. 
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